The J-triplet Cooper pairing with magnetic dipolar interactions 
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Recently, cold atomic Fermi gases with the large magnetic dipolar interaction have been laser 
cooled down to quantum degeneracy. Different from electric-dipoles which are classic vectors, atomic 
magnetic dipoles are quantum-mechanical matrix operators proportional to the hyperfine-spin of 
atoms, thus provide rich opportunities to investigate exotic many-body physics. Furthermore, un- 
like anisotropic electric dipolar gases, unpolarized magnetic dipolar systems are isotropic under 
simultaneous spin-orbit rotation. These features give rise to a robust mechanism for a novel pairing 
symmetry: orbital p-wave (L = 1) spin triplet {S — 1) pairing with total angular momentum of the 
Cooper pair J — I. This pairing is markedly different from both the '^He-B phase in which J = 
and the '^He-A phase in which J is not conserved. It is also different from the p-wave pairing in the 
single-component electric dipolar systems in which the spin degree of freedom is frozen. 

PACS numbers: 03.75. Ss, 74.20.Rp, 67.30.H-, 05.30 
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Ultracold atomic and molecular systems with electric 
and magnetic dipolar interactions have become the re- 
search focus in cold atom physics [l|-[3| ■ When dipole mo- 
ments are aligned by external fields, dipolar interactions 
exhibit the d^^-^z^-type anisotropy. The anisotropic 
Bose- Einstein condensations of dipolar bosons (e.g. ^^Cr) 
have been observed [sl-KLlj. For the fermionic electric 
dipolar systems, ''°K-®''Rb has been cooled down to 
nearly quantum-degeneracy Q . Effects of the anisotropic 
electric dipolar interaction on the fermion many-body 
physics have been extensively investigated. In the Fermi 
liquid theory, both the single particle properties and col- 
lective excitations exhibit the dr'^-^z^ anisotropy [12l - ll7l |. 
In the single-component Fermi systems, the leading or- 
der Cooper pairing instability lies in the p-wave channel, 
which is the simplest one allowed by Pauli's exclusion 
principle. The anisotropy of the electric dipolar inter- 
action selects the instability in the p^-channel, which is 
slightly hybridized with other odd partial wave channels 
[2fl-[23|. For two-component cases, the dipolar interac- 
tion leads to anisotropic spin-triplet pairing, and its or- 
bital partial wave is again in the p^-channcl [28 31]. The 
triplet pairing competes with the singlet pairing in the 
hybridized s + dr^ -^z^-chsLunel. The mixing between the 
singlet and triplet pairings has a relative phase ±-|, which 
leads to a novel time-reversal symmetry breaking Cooper 
pairing state [29l |. 

An important recent experimental progress is the laser 
cooling and trapping of magnetic dipolar fermions of 
^^^Dy and ^^'^Dy with large atomic magnetic moments 
(lOfis) [il, H- There are important differences be- 
tween magnetic and electric dipolar interactions. Elec- 
tric dipole moments are essentially non-quantized clas- 
sic vectors from the mixing between different rotational 
eigenstates, which are induced by external electric fields 
[Sj .4], thus electric dipoles are frozen. In the absence 
of external fields, even though at each instant of time 
there is a dipole moment of the heteronuclear molecule, 
it is averaged to zero at a long time scale. In contrast, 
magnetic dipole moments of atoms are intrinsic, propor- 
tional to their hyper-fine spins with a Lande factor. Un- 



polarized magnetic dipolar Fermi systems are available, 
in which dipoles are defrozen as non-commutative quan- 
tum mechanical operators, thus lead to richer quantum 
spin physics of dipolar interactions. Furthermore, the 
magnetic dipolar interaction is actually isotropic in the 
unpolarized systems. It is invariant under simultaneous 
spin-orbit rotations but not separate spin or orbit rota- 
tions. This spin-orbit coupling is diiferent from usual 
single particle one, but an interaction effect. It plays an 
important role in the Fermi liquid properties such as the 
unconventional magnetic states and ferro-nematic states 
predicted by Fregoso et al [l^ [l^ . 

It is natural to expect that magnetic dipolar interac- 
tion brings novel pairing symmetries not studied in con- 
densed matter systems before. The systems of ^^^Dy and 
^^•^Dy are with a very large hyperfine spin of F = ^, 
thus their Cooper pairing problem is expected to be very 
challenging. As a first step, we study the simplest case 
of spin- 2, and find that the magnetic dipolar interac- 
tion provides a novel and robust mechanism to the p- 
wave (L — 1) spin triplet {S — 1) Cooper pairing to 
the first order of interaction strength, which comes from 
the attractive part of the magnetic dipolar interaction. 
In comparison, the p-wave triplet pai ring in usual con- 
densed matter systems, such as ^He |32h34| . is due to 
the spin-fluctuation mechanism, which is at the second 
order of interaction strength (see Refs. [33, |3g| for re- 
views) . This mechanism is based on strong ferromagnetic 
tendency from the repulsive part of the '^He-'^He inter- 
actions. Furthermore, the p-wave triplet Cooper pair- 
ing symmetry patterns in magnetic dipolar systems are 
novel, which do not appear in '^He. The orbital and spin 
angular momenta of the Cooper pair are entangled into 
the total angular momentum J = 1, which is denoted as 
the J-triplet channel below. In contrast, in the ^He-i? 
phase [33|, L and S are combined into J — 0; and in 
the ^He-A phase, L and S are decoupled and J is not 
well-defined [3 21 |3J] . There are two competing pairing 
possibilities in this J-triplet channel with different val- 
ues of Jz'. the helical polar state {Jz = 0) preserving time 
reversal (TR) symmetry, and the axial state {Jz = ±1) 



breaking TR symmetry. The helical polar state has point 
nodes and gapless Dirac spectra, which is a time-reversal 
invariant generalization of the '^He-A phase with entan- 
gled spin and orbital degrees of freedom. In addition to 
usual phonon modes, its Goldstone modes contain the 
total angular momentum wave as entangled spin-orbital 
modes. 

We begin with the magnetic dipolar interaction be- 
tween spin-i fermions 

Vai3,P'a'{r) = —{Saa' ' Spp' - 3{Saa' ' r){Spi3' • r) } , (1) 

where r is the relative displacement vector between two 
fermions; /i is the magnitude of the magnetic moment. 
Such an interaction is invariant under the combined 
SU{2) spin rotation and 5*0(3) space rotation. In other 
words, orbital angular momentum L and spin S are not 
separately conserved, but the total angular momentum 
J = L + S remains conserved. Its Fourier transformation 
reads ^ 

47r -> -^ -> -> 

VQ/3;/3'Q'(q) = ^M^ {3(5*00' • q){Sl3P' ■ q) - Saa' ■ Spi3'}.{2) 

The Hamiltonian in the second quantization form is writ- 
ten as 



H = Y,Hk)-fIc]cUk)Ca{k) + ^ 



2V 



k.a 



J2 Vap;p'a'ik - k')Plp{k;q)Pp'c.'ik';q), (3) 

k,k' ,q 

where e(fc) — hk"^ /{2'm); fXc is the chemical potential; 
Pp'a'{k]q) = ci3'{—k+q}ca'{k+q) is the pairing operator; 
the Greek indices a, /3, a' and /3' refer to f and |; F is 
the volume of the system. We define a dimensionless 
parameter characterizing the interaction strength as the 
ratio between the characteristic interaction energy and 

the Fermi energy: A = Eint/Ep = f ^^2*^2^ ■ 

We next study the symmetry of the Cooper pairing in 
the presence of Fermi surface, i.e., in the weak coupling 
theory. An important feature of the magnetic dipolar in- 
teraction in Eq. ([T]) is that it vanishes in the total spin 
singlet channel. Thus, we only need to study the triplet 
pairing in odd orbital partial wave channels. Consid- 
ering uniform pairing states at the mean-field level, we 
set g = in Eq. ([3]) , and define triplet pairing opera- 
tors Ps{k), which are eigen-operators of Siz + S2Z with 
eigenvalues Sz — 0,±1, respectively. More explicitly, 
they are Po{k) = j^[Pn{k) + Pitik)], Pi{k) = Pn{k), 

P-i{k) = Pii{k). The pairing interaction of Eq. ([3]) 
reduces to 



Hr, 



1 

2V 



J2 {KT.i(fc;fc')^l(fc)^<(fc')},(4) 



where 

aPl3'a' 
{Va/3.,P'a'ik-k')-Vap,P'a'ik + k')}. (5) 



{lsz\^0i^/3) is the Clebsch-Gordan coefficient for two 

spin-i states to form the spin triplet; and Vs^s'^{k;k') 

is an odd function of both k and k' . 

The decoupled mean-field Hamiltonian reads 



H, 



mf 



-E 

2V ^^ 

k 



- 7^> '*Hfc) 



e(fc)/ ^ap{k) 

mi 



A^Jfc) 



*(^),(6) 



where we only sum over half of the momentum space; 
^{k) = e{k) — ^ch and ^c/i is the chemical potential; 
*(fc) = (c^(fc),c;(fc),4(-fc),c];(-fc))^; A„^ is defined 

as Aq;3 = '^sA^^A\'^\P)* ^s^- As^ satisfies the mean- 
field gap function as 



A-(^^) = ^E 

k' .s'^ 



Vl,,^{k;k'){\P^,Sk')\) 



VlAk;k')[K{k')~—]/\,,{k'), (7) 
(/TTJ'^ ' Zek 

where K{k') = ta,nh[§E,{k')]/[2E,{k')]. The integral in 
Eq. d?]) is already normalized following the standard pro- 
cedure [20]. For simplicity, we use the Born approxima- 
tion in Eq. ([7|) by employing the bare interaction poten- 
tial rather than the fully renormalized T-matrix, which 
applies in the dilute limit of weak interactions. The pair- 
ing symmetry, on which we are interested below, does 
not depend on the details that how the integral of Eq. 
^ is regularized in momentum space. The Bogoliubov 



quasiparticle spectra become Ei^2{k) = \/S.k + ^i.2ik), 
where A^ 2{k) are the eigenvalues of the positive-definite 

Hermitian matrix A^{k)A{k). The free energy can be 
calculated as 

F = -I y In [2 cosh ^1 

k,i=l,2 

\- J2 {XAk)V^Xik;k')^,,{k)}A&) 



2V 






rT-l, 



k,k'. 



where Vg ^, (fc; fc') is the inverse of the interaction matrix 
defined as 

^ E vls,ik;k')v^;^(k';k")^S^j:,A^,.,. (9) 



We next linearize Eq. ^ around Tc and perform the 
partial wave analysis to determine the dominant pairing 
channel. Since the total angular momentum is conserved. 



we can use J to classiiy the eigen-gap functions denoted 
as t/)""'^'^' (fc). The index a is used to distinguish different 
channels sharing the same value of J. (/>^''^'^^(fc) satisfies 

dflk 



No 



An 



K^,, {k; k')r/'' {k') = w'^jct^f'^ (fc), (10) 



™ -^ is the density of state at the Fermi 



where Nq - ^^5^ 

surface; w" are dimensionless eigenvalues; fc, k' are at the 
Fermi surface. Then Eq. ([7|) is linearized into a set of 
decoupled equations 



0a;jj.|^ + w5[ln(2e^w)/(^fcBT)]} = 0, 



(11) 



where Co is an energy scale at the order of the Fermi 
energy playing the role of energy cut-off from the Fermi 
surface. 

The decomposition of V^^^, (fc; k') into spherical har- 
monics can be formulated as 



Lm.L'in' 



VLms^;L'.n's'Yl^{nk)YL>ra'{^T:,), (12) 



where L = L' oi L — L' ± 2, and L,L' are odd num- 
bers. The expressions of the dimensionless matrix el- 
ements VLms^-L'm's'^ are lengthy and will be presented 
elsewhere. By diagonalizing this matrix, we find that the 
most negative eigenvalue is w"'^^ = — 37rA/4 lying in the 
channel with J — L — 1. All other negative eigenval- 
ues are significantly smaller. Therefore, dominate pair- 
ing symmetry is identified as the J-triplet channel with 
L = S' = 1 in the weak coupling theory. Following the 
standard method in Ref. [2(|, the transition tempera- 

2e''Ld-TZ7 



ture Tc is expressed as Tc 



For a rough 



estimation of the order of magnitude of Tc, we set the 
prefactor in the expression of Tc as Ef . 

In order to understand why the J-triplet channel is 
selected by the magnetic dipolar interaction, we present 
a heuristic picture based on a two-body pairing problem 
in real space. Dipolar interaction has a characteristic 
length scale adp = m^"^ /It? at which the kinetic energy 
scale equals the interaction energy scale. We are not 
interested in solving the radial equation but focus on the 
symmetry properties of the angular solution, thus, the 
distance between two spins is taken fixed at adp. We 
consider the lowest partial-wave, p-wave, channel with 
L = 1. The 3x3 = 9 states {L = S = \) are classified 
into three sectors of J = 0, 1 and 2. In each channel of 
J, the interaction energies are diagonalized as 



En ~ E, 



dp. 



1 

E, - --Ed 



1 

E2 - Y^Edp, 



(13) 



respectively, where Edp = /J' /cidp- Only the total angu- 
lar momentum triplet sector with J — 1 supports bound 
states, thus is the dominant pairing channel and is con- 
sistent with the pairing symmetry in the weak-coupling 
theory. 




FIG. 1: The spin configurations of the two-body states with 
a) J = 1 and jz = and b) J = jz — 0. The interactions are 
attractive in a) but repulsive in b). 



This two-body picture applies in the strong coupling 
limit. Although a complete study of the strong coupling 
problem is beyond the scope of this paper, this result pro- 
vides an intuitive picture to understand pairing symme- 
try in the J-triplet sector from spin configurations. We 
define that x^ andp^(f2) are eigenstatcs with eigenvalues 
zero for operators e^ • {Si + S2) and e^ • L (/i = x, y, z), 
which are the total spin and orbital angular momenta 
projected along the Cu-direction. The J-triplet sector 



states are 0^(rj) 



V2 l^' 



uxXvPxi^) with (/)^ satisfying 



(e^ • J)4)p = 0. For example. 



0,(0) 



1 

XCOSt. 



XxPyi^) -XyPxi^)] 

sin0{|aeji|a,,)2 + |%)i|%)2},(14) 
+ ysm(j> and \ae^) and |/3ep) are eigen- 



where 1 

states oicp-a with eigenvalues of ±1. As depicted in Fig. 
[1] (a), along the equator where (/), has the largest weight, 
two spins are parallel and along f, thus the interaction is 
dominated by attraction. On the other hand, the eigen- 
state of J = reads 

M^) = X^.P^.m = ■^{\an)iWn)2 + |A2)i|«o)2},(15) 

where |aa) and \/3q) are eigenstatcs oi (l ■ a with eigen- 
values ±1. As shown in Fig. [T](b), along any direction of 
f2, two spins are anti-parallel and longitudinal, thus the 
interaction is repulsive. 

Let us come back to momentum space and study the 
competition between three paring branches in the J- 
triplet channel under the Ginzburg-Landau (GL) frame- 
work. We define 

A,(fc) = -^[-Ai{k) + A-i{k)], 
v2 



Ay{k) 



^/2 
A,(fc) = Ao(fc) 



Ai(fc) + A_i(fc)], 



(16) 




FIG. 2: The angular distribution of the gap function |A(fc)|^ 
v.s. cos 9k in the helical polar pairing state (the red line) and 
the axial pairing state (the black line). 



The bulk pairing order parameters are defined as A^ — 

V Sfc ^/i^M(^)' 'where no summation over /i is assumed. 
We define pairing parameters and their real and imagi- 
nary parts as the following 3-vectors A = (A^;, Aj,, A^). 
The GL free energy is constructed to maintain the C/(l) 
and 5*0(3) rotational symmetry as 



F 



aA* -A + 71 1 A* 



where 



NoHy)- 



72|A*xA|2, (17) 



(18) 



The sign of 72 determines two different pairing structures: 
ReA II ImA at 72 > 0, and ReA _L ImA at 72 < 0, re- 
spectively. Using the analogy of the spinor condensation 
of spin- 1 bosons, the former is the polar pairing state and 
the latter is the axial pairing state |37h40| . 

For the polar pairing state, the order parameter con- 
figuration can be conveniently denoted as A = e*"^|A|z 
up to a U{1) phase and S'0(3)-rotation. This pairing 
carries the quantum number J^ = 0. The pairing matrix 
A^'^ = i|A|[fcyCri - k^cr2)i(72]ai3 reads 



A^ 



-|A| 



\ y ' ^'^x) 








(19) 



It equivalents to a superposition of Px T iPy orbital con- 
figurations for spin-ft (ii) pairs, respectively. Thus, 
this pairing state is helical. It is a unitary pair- 
ing state because A^^A is proportional to a 2 x 2 
identity matrix. The Bogoliubov quasiparticle spectra 
are degenerate for two different spin configurations as 



^k a ~ \/ik~^ |AP'(fc)p with the anisotropic gap func- 
tion |AP'(fc)|2 = ilApsin^^fc depicted in Fig. [^ They 
exhibit Dirac cones at north and south poles with oppo- 
site chiralities for two spin configurations. 

Similarly, the order parameter configuration in the ax- 
ial pairing state can be chosen as A = -y=e^'^\A\{ex+iey) 
up to the symmetry transformation. This state carries 
the quantum number oi Jz = 1. The pairing matrix 




1.12 



1.08 



1.04 



1.00 



FIG. 3: The ratio of the angular integrals of the free en- 
ergy kernels 2/(gm, -flm)' which is always larger than 1. This 
means that the polar pairing is favored at the mean-field level. 



a/3 - i^\^\{[kz{cfi + ia2) + (Jz(kx 



A ax 

^Q/3 — 

the form 



+ iky]ia2}ap takes 



Ar^ = ^|A| 



2\kx ^ky) 



2 ("'K + "^ky) 







(20) 



This is a non- unitary pairing state since A^^A = 
|A|2[i(l + kl) + kz{k ■ a)]. The Bogoliubov quasi- 
particle spectra have two non-degenerate branches 
with anisotropic dispersion relations as E1\(k) — 



^i + |Ar 



lA^f 



(fc)p. The angular gap distribution 

Ap(l ± cosSfc)^ is depicted in Fig. 



[21 Each of branch 1 and 2 exhibits one node at 
north pole and south pole, respectively. Around the 
nodal region, the dispersion simplifies into i?i,2(fc) = 
\ v'l{kz T kfY + ^|Ap(fc||/fc/)'', which is quadratic in 

the transverse momentum k\\ = \/kx + k^. 

At the mean-field level, the helical polar pairing state 
is more stable than the axial state. Actually, this con- 
clusion is not so obvious as in the case of ^He-B phase, 
where the isotropic gap function is the most stable among 
all the possible gap functions [33|. Here, the gap func- 
tions are anisotropic in both the polar and helical pairing 
phases. We need to compare them by calculating their 
free energies in Eq. ([S]) . The second term contributes the 
same to both pairing phases. Thus, the first term deter- 
mines the difference in free energies. Let us define the 
ratio between angular integrals of the free energy kernels 
in Eq. ([8|) of the two phases as 



2/(Ai,A2) 



/ dQk 2 In 



2cosh|./er 



|Ap'(fc)|^ 



/df^fcl]±ln 



2cosh^A/C^ 



|A|^(fc)|2 



,(21) 



where Ai = 



2/(Ai,A2) is numerically 



1 \ _ 1 

plotted in Fig. [3l For arbitrary values of /3, £_k, and |A|, 
y is always larger than 1. Therefore, the polar state is fa- 
vored more than the axial state. This can be understood 



from the convexity of the nonlinear term in Eq. ([8]), 
which favors isotropic angular distributions of |A(fc)p 
[4a | . Although neither gap function of these two states is 
absolutely isotropic as in the ^He-B phase, the polar gap 
function is more isotropic from Fig. [2]and thus is favored. 
However, we need to bear in mind that we cannot rule 
out the possibility that certain strong coupling effects can 
stabilize the axial state. In fact, the '^He-A phase can be 
stabilized under the spin feedback mechanism [35| , which 
is a higher order effect in terms of interaction strength. 

Next we discuss the classification of Goldstone modes 
and vortices in these two states. In the helical polar 
state, the remaining symmetries are SOj{2) x Z2 as well 
as parity and time- reversal (TR), where Z2 means the 
combined operation of rotation n around any axis in the 
xy-plane and a flip of the pairing phase by n. The Gold- 
stone manifold is 

[S0j{3) X Uc{l)]/[S0j{2) ® Z2] = [Sj X Uc{l)]/Z2(22) 

The Goldstone modes include the phase phonon mode 
and two branches of spin-orbital modes. Vortices in this 
phase can be classified into the usual integer vortices 
in the phase sector and half-quantum vortices combined 
with TT-disclination of the orientation of A. In the axial 
state, the rotation around z-axis generates a shift of the 
pairing phase, which can be canceled by a C/c(l) trans- 
formation, thus, the remaining symmetry is SOj^-^{2). 
The Goldstone manifold is S'^ x Uc{l). Only integer vor- 
tices exist. 

In summary, we have found that the magnetic dipolar 
interaction provides a robust mechanism at first order in 
the interaction strength for a novel p-wave (L = 1) spin 



triplet {S = 1) Cooper pairing state, in which the total 
angular momentum of the Cooper pair is J = 1. This 
is a novel pairing pattern which does not appear in ^He, 
and, to our knowledge, neither in any other condensed 
matter systems. These pairing states include the TR 
invariant helical polar pairing state and the TR breaking 
axial pairing state, both of which are distinct from the 
familiar ^He-A and B phases. 

Many interesting questions are open for further explo- 
ration, including the topological properties of these pair- 
ing states, vortices, spin textures, and spectra of collec- 
tive excitations. The above theory only applies for spin- 2 
systems, in which the magnetic dipolar interaction is too 
small. For the pairing symmetry in a magnetic dipolar 
system with a large spin S, our preliminary results show 
that the basic features of the J-triplet pairing remains. 
The spins of two fermions are parallel forming Stot ~ 2S 
with orbital partial- wave L = 1, and the total J = 25*. In 
the current experiments in Ref. [4l|, the highest attain- 
able density reaches 4 x lO^^cm"-"^ for -"^^-"^Dy atoms with 
S = ^. The corresponding dipolar energy is Eint ~ 2nK 
and the Fermi energy for unpolarized gases Ef w 13.6nK, 
and thus A = Eint/Ef sa 0.15. If we use the same for- 
mula of w^^^ above for an estimation of the most nega- 
tive eigenvalue, we arrive at Tc/Tf w 0.06, which means 
that Tc « 0.8nK. Although it is stiU slightly below the 
lower limit of the accessible temperature in current ex- 
periments, we expect that further increase of fermions 
density, say, in optical lattices will greatly increase Tc- 
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